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Abstract. Radiation damage induced by ion beams is traditionally treated at different levels of theoretical approaches, for the 
different scales and mechanisms involved. We present here details of a combined approach that, from a method at a nanoscopic 
scale, attempts to merge with higher scales existing results, by "tuning" the analytical method employed when extended to 
larger scale and so yielding a consistent picture of the entire process. Results will show the possibility to get a good agreement 
with macroscale methods and, on the other hand, to produce a reliable electron energy spectra to be used for DNA damage 
estimations. 
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>»: INTRODUCTION 

Presently, ion-beam cancer therapy, after a number of outstanding clinical results, is a widely used treatment tool 
based on modern theoretical and computational methods of basic science fl]. Indeed, while the general principles 
governing the phenomenon are well understood for decades, many details of the involved physical processes, that lead 
from the incidence of an energetic ion in the media to the biological damage, are still far from being explained and 
quantitatively estimated at a nanoscopic level. More, between the many different approaches exploited at different 
\^ , scales, with a number of very accurate works |2, 3, 4], there is no unified picture of the entire mechanism . A new 

Ps| ' approach is presented for a preliminary estimation of the damage occurring in a DNA unit placed at a given distance 

from a ^^C ion track (see also in this proceedings book |5]). A major part of the damaging effects following the 
f-^ ■ passage of an energetic ion in a tissue like medium, are connected, directly or indirectly, to the number and energy 

00 ! distributions of secondary electrons produced by the primary ionization of the medium molecules . Our method starts 

^^ ■ from a calculation of the energy spectrum of the electrons originally produced by the ion beam [3,4], followed by the 

evaluation of their possible pathways from the track. This energy spectrum is obtained from the analytical calculation 
of the singly differentiated cross section of ionization of the medium. After integration this yields the penetration 
depth and the typical profile of the Linear Energy Transfer (LET). The key idea is then to set up the equation for 
electron production , and then integrate it to large scale for checking consistency, with macroscopic scale results. In 
other words, by looking at the macroscopic quantities one can adjust the model for electron production the analytical 
method in order to produce consistent results from the nanoscale to the millimetric scale. To the best of our knowledge 
there are only a few attempts calculating dose profiles in a similar way |6] Possible extensions of the model that are in 
progress will be also mentioned. 

THE MACROSCOPIC SCALE: ENERGY DEPOSITION OF THE ION IN THE 

MEDIUM 

The key quantity of interest for the ionization process, linking the electron production and the propagation of the ion 
in the medium, is da/dW, the Singly Differentiated Cross Sections (SDCS). As we showed in our previous papers 
|l3,|8[l, it is possible to obtain it by a semiempirical formula by Rudd |9[], similar to the Bethe-Bloch formula for energy 
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FIGURE 1. Total Cross section of ionization . In the inset the Zeff behaviour as a function of T according to llOll 



loss, with the inclusion of corrections for imposing a correct asymptotic behaviour.; 



d(y{W,T) 
dW 



Ni 



'-rp 



m T 



(1) 



Where / is a generic function, with a number of fitted parameters, of W, the kinetic energy of emitted electrons and 
T, the incoming projectile energy per nucleon (E /u), scaled by /, the ionization potential of a given shell /, with 
occupation number A^,, and where m/M is the ratio between the electron and ion masses, and z the ion's charge(for the 
details seefS]). 

The correct value of z, different from the initial bare ion charge, is one of the most important correction to include. 
One has indeed , to take into account the physical effect of electron capture, experienced by the ion in the last part of 
its trajectory, when its velocity is relatively small, that causes the initially bare ion to collect a number of electrons and 
to reduce its effective charge. For the total ionization cross section, obtained by integration on the secondary electron 
energies 

•d(7(W,T) 
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70 dW 

and plotted in Fig. 1, we can see that this effective charge, computed for an ion at velocity j5c according by a formula 
reported by Barkas lUOl] (inset of fig. 1) 

z,//=z(l-exp(-125i3z-2/3)), 



(3) 



undergoes a drastic decrease till about one half of its initial value, just in the region where the cross section is maximum 
(Bragg peak), that is the most important for damaging studies. This suggests experiments that could mimic the Bragg 
peak scenario by easiest conditions, that is by simply using a half-stripped carbon (C+^) and a collision energy 
approaching 100 keV. Another interesting effect related to this charge variation, consists in a displacement of the 
maximum of the curve from 0. 1 to 0.3 MeV/u, as can be seen in fig. 1 . 

But the most stringent test of the used formula is, after the integration, to retrieve the Linear Energy Transfer. An 
energy loss by the ion in a single collision, W + U , integrated with the SDCS and multiplied by the water density n 
gives the Linear Energy Transfer (LET, dT/dx) as a a function of the penetration depth x that shows, on a millimetric 
scale,the typical Bragg peak profile: 
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FIGURE 2. Linear Energy Transfer profiles near the Bragg peak for C''^ 
data(*) are shifted by 7.2 mm,in order to compare the shapes 



at 400 MeV/u obtained with different methods. Our 



where Tq is the initial ion energy per nucleon of incidence of the ion. That allows us to compare (fig. 2) with 
experimental data [Ulll and simulations with Monte Carlo model for Heavy Ion Therapy (MCHIT) based on the 
GEANT4 toolkit Ill2ll . We have already shown in our previous paper| 8 ] the importance of the inclusion of a relativistic 
correction, that shifts the Bragg peak position up to the 40% of its value. After that the longitudinal deviation in the 
depth of the Bragg peak is within the 3% of the whole penetration depth, a value that we can argue to be due to the 
energy spent for excitation without ionization by the particle. The accounting of this quantity is presently in progress. 
The introductionof the effective charge z^jf doesn't affect the position of the peak but reduces strongly its height as 
seen for the C7(r). Anyway this reduction is still not sufficient for yelding the correct height and shape of the peak as 
shown in fig. 2 (thin curve), where, for better clearity, the remaining displacement in penetration depth is compensated 
by properly shifting our curve. To this end our calculation does not include the stochastic nature of ion energy loss 
which appears both in experiments and MCHIT results. Then we introduce such stochasticity by accounting for the 
energy straggling due to multiple Coulomb scattering of ion in water 1 13]. A gaussian spreading of the LET according 
to a longitudinal-straggling standard deviation parameter Ox — Q-^inm computed by HoUmark et al. Iil4il . for carbon 
ions of therapeutic energies, gives an average value 
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where xq is a maximum penetration depth of the projectile. With this correction we got perfect agreement with a 
Monte Carlo simulation where the nuclear fragmentation channels are switched off, that verifies our model in the 
absence of fragmentation (dashed curves in fig. 2). We can then simulate the nuclear fragmentation by accounting the 
fluence attenuation of primary ions due to nuclear reactions along the pathway. The total nuclear charge changing 
cross section is defined for nuclear fragmentation processes when the initial nucleus of charge Z=6, changes to any 
other nucleus with Z<6. This value is measurable experimentally and we take the value reported by Chu for incoming 
energies range consistent with our case) f 15] : CTz<6 ~ 1248 mb and connected to the mean free path 



^nuc(To) — 
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where the mean free path is slowly dependent on the initial energy per nucleon To, and could be assumed constant in 
the therapeutical range (from 200 to 450 MeV/u). Thus we can consider the attenuation of the fluence <I>c as a factor 



k[x) = <I>c(x)/4>c(0) = exp{-x/l„,,c). 



(8) 
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FIGURE 3. Volumetric energy deposition along and around the track for C^^ at 400 MeV/u. Each countour step stays for an 
increase of 1 MeV/mm^. 

By including this factor to previous to eq|6] we obtain the effective energy transfer due to primary ions. As evident in 
figure 2(solid curve), the assumption not to take into account LET components from the nuclear secondary fragments 
will lead to underestimation in the region before and after the Bragg peak, while the estimation of the Bragg peak 
height, where the major contribution comes from the specific C ion, is consistent with experimental data and MCHIT 
results. 

While the analytical method is computationally very fast, our outcome could be used to easily compute also the 
volumetric energy deposition in the medium, by accounting not only the energy straggling in the propagation direction, 
but also in the radial components by an empirical formula provided by Chu depending on the charge and mass of the 
projectile llSll . 



a,{x;Z,A) = 



Z'A'i 



(9) 



where the empirical parameters are respectively a=0. 00373, b=0.896, c=0.207, d=0.396. 

It is then possible to map the energy deposition as in figure 3. In that contour plot we can easily see that the scale of 
this lateral spread due to stochastic effects is of the order of several millimeters. 

The experimental procedure usually accounts for a cylindrical spread of the initial beam with a Full Width at 
Half Maximum (FWHM) of the order of 3 mm, that is here also included by an additional lateral spreading term 
(do = FWHM/{2V2ha) 
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where S is the unit of surface for dimensional consistency. 



THE NANOSCOPIC SCALE: SECONDARY ELECTRONS AND DNA DAMAGE 



Once tested the reliability of the method at large scale, we can extract from the cross sections results information on 
the damaging electrons that are produced. First of all we want to investigate whether the electron showers produced 
by different single ions have negligible overlap between each other. As it is known from Monte Carlo simulations il6il 
and experimental data lUTIl the penetration lengths of the electrons of that energy range(I-100 eV, as shown below) in 
water are within 10 nm, we can take as the effective radius of the track, r„, this value. Thus considering a typical ion 
flux used as / = lO^ions cm^^s^' and the whole irradiation time, usually 20 min, we get the average number of ions 



passing an area of the size of the effective track as 

N = nrllt (11) 

and assuming a stochastic process, governed by a Poisson statistics we get the probability of k hits in the same area as 

P{k}^e^^ (12) 

that provides already for k=2, P(2)=5x 10^^, when considering the total irradiation time and P(2)= 1.8 x 10^^ when 
taking t=l min, time sufficient for repairing mechanisms to start and then levelling out the previous consequences of 
the damage. We can then consider the damaging effects connected to different tracks to be independent. 

One can estimate the damage produced by ion tracks by considering values plotted in figure 1. In the Bragg 
peak region Gbp ~ (7(0.3) w 40A^, then the number of electrons produced for 1 nm of the ion's trajectory will be 
(multiplying by the density of water n — 33A nm^^) equal to naBP ~ 1 3 nm^ ^ . This quantity divided by the effective 
area along the track, as mentioned above, gives the number density of the electron plasma along the track: 

Me = — ^ w 0.(H5nm-^ (13) 

As it is known the conventional measure of the DNA damage is related to the number of double strand breaks, that 
is the occurrence of 2 damaging events on opposite strands within a single DNA convolution of length lDNActmv=^-^ 
nm and radius roNA=^-^ nm. Now, considering the volume of that convolution VoNAamv = ^foNA^DNAconv ~ 12.9 
nm^ we get the average number of electrons in such a volume, within the effective volume around the track, as 
f^e/DNA ~ ^DNAcom>ne = 0.6. This number reveals that the electron density along the track is on average remarkably 
high compared to the minimum density requested to have a DSB. We can replace this average number by a distribution 
of the radial density of this electron plasma as a function of the distance from the track, still expressed in units of 
volume of a DNA convolution, by imposing a gaussian attenuation from the track as 

rie/DNAir) = ^exp-{p- /2al) (14) 

with Oa = '"a/3 and A'^, = raU^iDf^j^ and we can see from this profile (figure 4) that at short distances the density is close 
to 2. This very high density, is a peculiar effect of carbon ions, while e.g. for protons, as mentioned above, the much 
smaller cross sections will give much lower density. So in this case the simultaneous damage due to different electrons 
on the same convolution, will be emphasized, while lower densities could produce a double strand break only by the 
action of the same single electron. That is another reason why we are going to focus our interest on this special kind 
of damage. Heavier ions, of course, with their larger charge would produce even larger densities, but their use shows 
different contraindications yj], in particular related to increased nuclear fragmentation. 

For quantitative estimations, one has to take into account the different trajectories of the produced electrons. From 
the SDCS of eq.O we can extract the spectrum of the energies of the produced secondary electrons in a given A;ic 
segment 

dW dW ^ ' 

We show in figure 5 (left panel) these data for Ax =1 nm , and one can clearly see that, at any ion energy considered, 
the number of produced electrons with W>100 eV is negligible, and a large majority is produced below 30eV. This 
result is confirming the importance of the studies of biological effects induced by low energy electrons that are recently 
in progress both by experimental [ ISlli^l and theoretical 1120112 111 means. 

In the right panel of figure 5 it is shown the average energy of produced electrons according to the formula 

The behaviour of this quantity is characterized by a plateau around the Bragg peak, while at high T it reach a value 
around 60 eV. The results obtained with different method for the same process but considering liquid water Jlm are 




4 6 8 10 

distance from the track, r (nm) 

FIGURE 4. Number density distribution of the electron produced by a carbon ion track for a unit volume equal to a single DNA 
convolution. On the left, a scheme of the area around the track wandered by the secondary electron, and the DNA convoltion taken 
as unit. 
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FIGURE 5. (left panel) Energy distribution of secondary electrons produced by a C^+ ion, in a track length of 1 nm, for different 
values of T (MeV/u). (right panel) Average kinetic energy of produced electrons as a function of the ion energy compared to a 
similar calculation for liquid water Ii22i1 



also shown for comparison. Even if some deviations are evident, a qualitative agreement in the overall behaviour is 
achieved. 

The obtained electron distribution, dN/dW, can be profitably used for building models for quantitative estimation of 
DNA double strand breaks induced directly by these electrons or indirectly by radicals creation [5]. 

This approach shows the possibility, even with the use of many simplifications, to face the problem of the damage 
induced by ion beams in a quantitative way, by including, one by one, all the relevant physical effects appearing at 
different scales of space and time. We can point out then that there is plenty of room for experiments that may improve 
the understanding of any single part of these physical processes. 
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